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Abstract 

Circumbinary discs are commonly observed around post-asymptotic giant branch (AGB) systems and are known to 
play an important role in their evolution. Several studies have pointed out that a circumbinary disc interacts through 
resonances with the central binary and leads to angular momentum transfer from the central binary orbit to the disc. 
This interaction may be responsible for a substantial increase in the binary eccentricity. We investigate whether this 
disc eccentricity-puiuping mechanism can be responsible for the high eccentricities commonly found in post- AGB binary 
systems. 
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1. Introduction 

Post-AGB stars just left the AGB phase and evolve rapidly 
to hotter effective temperatures at constant luminosity, but 
they are not yet hot enough to ionise circumstellar material 
ejected during the AGB phase. Some post-AGB stars also 
show a near-infrared emission due to thermal radiation of 
hot circumstellar dust. The colour temperature of the dust 
is a good indication that there must by circumstellar dust 
close to the star and this SEP cha racteristic is always re- 
lated to binarity ()Van Winckel et a l. 2009, and references 
therein). The extr emely narrow CO emission lines (a few 



km s ^ width, e.g. iJura et ah 1995), th e presence of large 



(sub- micron) grains ( Gielen et al.ll2011 | ) and the detec tion 
of Oxygen-rich crystalline silicates ( Waters et al .111998 are 
distinct characteristics best explained by the long-term pro- 
cessing of dust in a stable circumbinary (CB) disc. 

However, a key problem still remains regarding the in- 
ability of population-synthesis models to reproduce the ob- 
served large eccentricit ies of post-AGB stars, similarly ob - 
served in barium stars (jPols et al.ll2003l : llzzard et al.ll2010t ). 
In the previous evolutionary phase (i.e. the AGB), the very 
large stellar radii lead to efficient circularisation of the 
binary orbit by tides ([Zahn 1977). Systems too close to 
accommodate an AGB star undergo Roche-lobe-overflow 
(RLOF) mass transfer, in most cases in an unstable regime 
which leads the system into a common envelope (CE). 
Due to friction between the stellar cores and the CE, the 
envelope may be ejected, and the orbit shrunk and cir- 
cularised. Post-AGB binaries with orbital periods shorter 
than 1, 000 d are expected to be circular whereas observa- 
tions reveal instead systems with periods from about 10^ d 
to 3, 000 d, often in eccentric orbits. A mechanism that in- 
creases the binary eccentricity is therefore required. 

If tides are as efficient as predicted (jZahnl |1977|) , a 
mechanism to increase the eccentricity must take place 
during the AGB or post-AGB evolution. Different mech- 
anisms have been investigated such as enhanced mass- 
transfer at periastron (i.e. differential mass loss dur- 



ing an orbital period; ISokeiJ |2000[ ). eccentricity pump- 
i ng induced by a wind- RLO F hybrid mass transfer 
(jBonacic Marinovic et al.| 200S ) and a kick to the newly- 
born white dwarf ( Izzard et al.ll201(il) . 

In this paper we investigate the possibility that CB discs 
cause the large eccentricities observed in some post-AGB 
systems. The current knowledge of CB discs is summarised 
in Section [21 In Section [3] we present our results and com- 
pare them to observed post-AGB periods and eccentricities. 
Section m discusses successes and potential problems of our 
model. Conclusions are drawn in Section [Sj 



2. Circumbinary discs 

In this section we enumerate the classes of stars known or 
suspected to be binaries that possess dust discs. 



2.1. Observed disc properties 

Dust discs are observed in very different classes 
of stars known to be binaries (jChesneaul 1201 if ). 
Stable discs a re observed around most known post- 
AGB binaries (IVan Winckell 120031 : Ide Ruvter eHI] 120061: 
IVan Winckel et aI]l2006^ . They are increasingly being iden- 
tified arou nd the central stars of planetary nebulae, e.g. 
NGC 2346 (jCosteroet al.lll986[ l. Other evidence exists such 
as eclipses from dust discs (iii NGC 2346, C PD-56°8032 
and M2-29. iHaiduk et al.ll2008t|Gesicki et al.| [20iQl. Some 
symbiot ic stars (lAngeloni et al.l^2007^ ■ hydrogen-deficient 
binary ("Netolick v et al.l 120091 ) and young stellar objects 
(Deroo et al. 2007al) are also known to possess long-lived 
CB discs. Around some silic ate J-type stars, the disc is 
resolv ed (iDeroo et al I20071J) . Some AGB stars (e.g. X 
Her ; iKahane fc Jural |1996[ ) have a very slowly expand- 
ing wind (a few km s~^), however none is a confirmed bi- 
nary system. Kepleria n discs are also observed in B[e] stars 
(jMeilland et al.ll2007T ): but only few B[e] stars are confirmed 
binaries and it remains unclear whether the B[e] behavior is 
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related to binarity (jLamers et alJlTOOSh . Fi nally, discs have 
been discovered around white dwarfs (e.g. iGansicke et al.l 
[2006) . Surveys for such discs around the o ldest white dwarfs 
have been unsuccessful (jKilic et al.ll2009f ). but around 15% 
of local white dwarfs show metal-rich material in their pho- 
tosp here, indicative o f accretion from a residual dust reser- 
voir (jSion et al.ll2009[ ) , although this reservoir may not nec- 
essarily originate from binary interaction. 

The variety of classes of objects that host CB discs 
and the wide range of orbital periods of the central stars 
(from 115.9 d for SA O 173329 - [va n Winckel et al. 200C| - 
to 2597 d for U Mon; iPollard fcCottrelll 11995) show that 
they are very common around binary systems and easily 
formed in late-type mass-transfer systems. The structure 
of discs surrounding youn g stellar objects and post- AGB 
stars is somewhat similar ()de Ruvter et al.ll2006r) . which is 
rather remarkable because of their (possibly) very different 
formation mechanism. 

The stable nature of the dust disc in post- AGB sys- 
tems can be inferred from the observation that none of 
them shows evidence for a current dusty mass-loss while 
the dust excess starts near dust-sublimation radius; the 
near-IR emitting material must be gravitationally bound in 
the system. In addition, post- AGB systems show depletion 
patterns, i.e. the gas in the disc is separated from the dust 
and s ubsequently reacc reted onto the star ([Waters et al.l 
119921 : iMaas et al.ll200a ). This re-accretion process occurs 
more efficiently if circumstellar gas and dust are trapped 
in a stable disc (fVan ^yinckel et all l2006f . Ide Ruvter et al.l 
l2006t[Gielen et aLll2008l ). Furthermore the presence of near- 
IR (from about 5 to 20 fim) excesses indicates that in 
all systems the circumstellar shell is not freely expanding 
but stored in a disc. The detection of cool (from about 
100 to 150 K) Oxyg en-rich crystalline silicate dust grains 
IjWaters et al.lll998r ) as well as the presence of large (sub- 
micron) grains are other indications of the long-lived natur e 
of these discs fe.g. lAngeloni et al.|[2007l: iGielen et al.ll2008[) . 
In some cases, Ke plerian rotation is detec ted at least in the 
inner disc region (jBuiarrabal et al.ll2005f ). 



2.2. Formation 

The formation of CB discs and their impact on the cen- 
tral star is not well understood. Most likely, CB discs in 
post-AGB stars remain from CE ejection or form during 
mass-transfer events in a binary system. The former case is 
supported by population synthesis models which show that 
most post-AGB systems that possess a CB disc are post CE 
systems. The possibility to form a CB disc from a fraction 
of the ejecte d CE has already been suggested from detailed 
CE models ( ISandguist et aP 119981: iNordhaus fc Blackmad 
120061: iKashr fc Soker 201 ij). It turns out that the CE is not 
necessarily completely ejected, so some gas remains around 
the system. While CE evolution is expected to lead systems 
to orbital periods of hours to hundreds of days, no disc has 
been observed around systems with periods shorter than a 
hundred days. 

On the other hand, many studies have investigated the 
possibility that matter flows thr ough the second Lagrangian 
point L2 during RLOF (e .g. Podsiadlowski et al.l 19921 : 
iFran kowski fc JorissenI l2007f . Ivan Rensbergen et al.l 1201 ih , 
possibly because of radiation pressure which could re- 
shape the Roche equipotentials by a reduc tion of the 
effective gravity of the mass-losing star (jSchuermanI 



I1972I : IFrankowski fc Tvlendal 1200 ll : iDermine et aH I20090 . 
Hydrodynamical simulations have shown that this can lead 
to th e formation of a CB disc in a close binary sys- 
tem (jSvtov et al.l l2009tl . However, the disc may not be 
stable and may disappear as soon as the mass transfer 
stops. Two-dimensional hydrodynamical simulations show 
that, for a typical slow and massive wind from an AGB 
star, the flow pattern is similar to Roche lobe over- 
flow, and that a small fraction of th e mass transferred 
to the companion flows through L2 (de Val-Borro et al.l 
120091 : iMohamed fc Podsiadlo^skil 120101) . The companion 
focuses the ejected material into the equatorial plane of 
the system (iTheuns fc JorissenI [l99l iTheuns et al.l 119961: 
iMastrodemos fc Morris! ri998|) . However, the possibility of 
forming a stable CB disc from a mass-losing giant star has 
not yet been confirmed by simulations. 



2.3. Evolution 

Surveys performed at wavelengths as long as 100 /im show 
that none of the post-AGB systems with a disc is observed 
to be surrounded by an outflow remaining from the AGB 
phase (JGielen et al.ll201lh . This means that the observed 
post-AGBs have left the AGB phase long enough ago for the 
last ejecta of the AGB star to reach a distance at which the 
emission at 100 /im cannot be detected. A lower limit of the 
disc lifetime can then be inferred. Using the code DUSTY 
(jlvezic fc Elitzuil[T997il for a post-AGB star with an effec- 
tive temperature of 4, 000 K, a luminosity of 5, 000 Lq and 
a wind composed only of silicate grains, we deduce a dis- 
tance of about 0.4 pc. Assuming a typical wind velocity of 
15 km s~^, this timescale is approximately 2.5 x 10'* yr. 

It is known that the CB-disc mass decreases during 
post-AGB and PN phases because material flows away from 
the o uter part of the dis c and reaccretes on the central 
stars. iGesicki et al.l (|2010|) pointed out the sharp decrease 
in CB-disc mass as the star ages. Finally, the discs observed 
around cooling white dwarfs have very low masses because 
all the gas was probably expelled and only the dust remains. 
An upper limit of the disc lifetime can then be estimated 
from the time since the mass transfer stopped, i.e. the end 
of the AGB phase, given by the timescale to cross the post- 
AGB and PN phases. While it has been est imated for sin - 
gle stars (of the order of 1, 000 - 10, 000 yr; lBlockeH[T995l ). 
this timescale is very uncertain for binaries. In single stars, 
the post-AGB phase ends as soon as the thin envelope re- 
maining from the AGB phase is entirely ejected. However if 
re-accretion occurs, e.g. from a CB disc, the envelope ejec- 
tion could take much longer, possibly up to 10^ yr. From 
the observed strongly depleted post-AGB stars in t he LMC 
(jRevniers fc van Wincke1l2007l : iGielen et al.ll2009f) . we can 
infer a dilution of their envelopes by a factor of about 10. As 
the typical envelope mass is about 10"'^ Mq (Blocker 1995), 
10~^ Mq must have been accreted from a CB disc. With a 
typical mass-loss rate of lO~^M0yr~^, the accreted mass 
therefore extends the post-AGB lifetime by about 10^ yr. 

The disc lifetime is det ermined by two compet ing mech- 
anisms: (i) disc accretion (jHartmann et al.lll998^ ■ and (ii) 
photoevaporation of the surface layer of the gas disc due 
to far- and extreme-ultraviolet radiation from the central 
star (jHoUenbach et al. 2000 ). For comparison, pre-main- 
sequence disc life times are thought to be 5 — 10 Myr 
(lYasui et al.ll2010D . 
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Once the star is on the white-dwarf cooling track, any 
furt her dissipation is set by the Poynting-Robertson drag 
(|Wei dcnschilling & Jackson 1993), with a time scale of the 
order of lO'' yr. The disc may remain detectable for about 
10® yr, long after any sign of the planetary nebula has dis- 
appeared. 



2.4. Resonant interactions 



iGoldrei ch fc Tremaind (|l979f ) and lArtvmowicz &: Lubowl 
(1994) describe the resonant and nonresonant interactions 
between a binary system and its CB disc using a lin- 
ear perturbation theory and assuming (i) the disc is thin 
(0.01 < H/R < 0.1, where H and R are the thickness and 
the half angular momentum radius of the disc respectively) 
and (ii) the nonaxisymmetric potential perturbations are 
small around the average binary potential. The binary po- 
tential is expanded in a series as 

$(r, 9, t) = I]™j0„,,K'^)exp [im(d- {l/m)Qbt)] , (1) 

of which only the real part is relevant, and I and m are 
integers, flh is the orbital angular frequency which, in an 

inertial frame, is fit, = ( ^3 ) , where a is the semi-major 
axis of the system, and M its total mass. The individual 
potential harmonics 4>m.i{i') rotate uniformly with pattern 
speed 

Vtp = {l/m)nb. (2) 

A resonance occurs in the disc when the pulsation k, of 
radial motion of a disc particle (on an epicyclic orbit with 
azimuthal pulsation fi) is commensurate with the angular 
frequency (fJ — Vip) in the frame rotating with component 
(Z, m) of the perturber. 



m{rt — Op) = ±K, 



(3) 



where the positive/negative sign corresponds to the 
outer/inner Lindblad resonance (LR) respectively. A coro- 
tation resonance (CR) occurs when Q. — Vtp. These resonant 
interactions result in the excitation of density waves located 
at radii 



r-CR = {m/lf/h 



and 



rLR = [(m±l)/Z] 



2/3, 



(4) 



(5) 



Many studies have investigated w hether these reso 
nances i ncrease the orbital eccentr i citv (lArtvmowicz et al. 



1991; Artvmowicz fc Lubowl 



2003). We apply the model of 



1994; Frankowski fc Jorissen 



Lubow fc Artvmowicd (|l996l ) 



for small and moderate eccentricities (e < 0.2). This model 
is based on results of SPH simulations which show that 
for a disc in which the kinematic viscosity, i/, is indepen- 
dent of the density, the torque is independent of resonance 
strength and width. The resonant torque can then be es- 
timated from the rate of change of the disc angular mo- 
mentum. The viscous torques within the disc balance the 
gravitational torque generated by resonances in the inner 
part of the disc and redistribute mass and angular momen- 
tum throughout the disc. The viscous torque is 



'h = Jd/r = Md^bV, 



(6) 



where the viscous timescale r = R^ /v — a ^{H/R) ^Vt ^ 
is of the order of 10^ yr, and Md, v and a are the mass, 



kinematic viscosity and the viscosity parameter of the disc 
respectively. The viscosity parameter is typically a ~ 0.1. 

The variation of the orbital separation due to the res- 
onant int eraction between the binary and the CB disc is 
given by (jLubow fc Artvmowicj[l996[ ) 

a Jd% 21 Md (hVcl 

— = I -p— = a — — ilh, (7) 

a Jorb"b m II \RJ R 

where Jorb is the orbital angular momentum and /i is the 
reduced mass of the binary system. The increase of the ec- 
centricity due to resonances can be given as a function of 
a/ a and depends on the binary eccentricity. At very small 
eccentricities (e < O.Iq^^^), the resonances weakly drive the 
eccentricity so the inner radius of the CB disc is maintained 
by the m = I resonances that take place close to the bi- 
nary at r/a < 1.7. The eccentricity pumping in that regime 
increases with e and is estimated (jLubow fc Artvmowicd 
Il996l ) to proceed as 



50e a 
a a' 



(8) 



up to a maximum of e = 5a ^'"^d/a at e ~ O.la^'^. At 
larger eccentricities (O.la^/^ ^ e ^ 0.2), the to = 2, / = 1 
resonance is the strongest contribution to the eccentricity 
driving and the disc is progressively pushed away as the 
binary eccentricity increases. In that regime, the eccentric- 
ity pumping decreases as 1/e. The eccentricity growth di- 
minishes as 1/e until e ~ 0.5 — 0.7, at which point reso- 
nances that damp the eccentricity begin to dominate (e.g. 
iRodig et al.ll201ir ). To include this effect, we assume the e- 
pumping to be inefficient (e — 0) for e > 0.7 (note however 
that a is not necessarily 0). The e— pumping efficiency for 
very small and moderate eccentricit ies (i.e. for e < 0.2) is 
shown in F ig. |T] and is given by (see iLubow fc Artvmowicz! 
I19961I2000D 



2 (1 - e-) 



I 
m 



-nfi;^ V"^ %/T 



(9) 



The inner radius of the disc, rin, is determined by the 
condition that the resonant torque, which tends to push 
the disc away, overco mes the viscous torque which acts 
to fill the inner region. lArtvmowicz fc Lubowl l|l994f ) esti- 
mate from SPH simulations the disc inner radius at differ- 
ent binary eccentricities and for different Reynolds num- 
ber of the disc gas, JZ — {H/R)^^a~^ (see Fig. 3 of 
lArtvmowicz fc Lubowl 119941 ) which we fit with 



nn(e,7^) = 1.7+-log(7^^/^)AU. 



(10) 



Finally, the variation of the orbital angular nionientuni 
in the centre-of-mass frame is given by 



'Jnrh — 'J< 



orb 



2a 



J. 



orb 



(i-e^y 



(11) 



Our model is in reason able (within a facto r two) agreement 
with the SPH results of lArtvmowicz et al.l (jl99ll ) , who give 
d/a ~ -4.3 X lO^'^nf,Md/Mb and ee ~ 1.9 x IQ-^nbMd/Mb 
at e == 0.1, for fi = 0.3, a = 0.1, 7^ = 10^, with a CB 
disc extending from r = 1.7a to r = 6a, with a surface 
density distribution a ex 7'^^. Our model gives d/a ~ — 2.8x 
IQ-'^nbMd/Mb and ee ~ 2.5 x IQ-^nbMd/Mb- 
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Figure 1. de/di in units of dlna/di as a function of e 
(Eq. E]) based on a model for very low and moderate eccen- 
tricities (jLubow fc Artvmowicdll996ll2000f ). for a circumbi- 
nary disc with a viscosity parameter a = 0.1. The curve is 
divided in three regions. For eccentricities from to 0.2, 
our model is strictly valid (red solid curve). For eccentric- 
ities from 0.2 to 0.7, the model is extrapolated outside its 
validity range according to an efficiency decreasing as 1/e 
(black dotted curve). At higher eccentricities, de/di = to 
account for the saturation behavior (blue dashed curve). 



2.5. Eccentricity gap 

Interestingly, pre-MS and MS systems lack nearly circu- 
lar orbits at long perio ds (called the eccentricity gap; 
see Mathieu, Il992l |1994[ ). whereas short-period systems 
are circularised by tides. Although the eccentricity dis- 
tribution of pre-MS and MS systems is likely the re- 
sult of numerous processes such as binary formation, stel- 
lar encounters and tides, interaction with a CB disc cer- 
tainly occurs as well. Discs commonly observed in pre-MS 
syste ms are similar to those observed around post-AGB 
stars (Ide Ruvter et al. _ 20061), wit h masses from 0.004 to 



0.3 M r7^ . lArtvmowicz et al.l lll99lh . 
' 19941) ■ iMathieu et all (|1995l) and 



Artvmowicz fc Lubow 



Lubow fc Artvmowicz 



19961) suggest that they pump the binary eccentricity. 



which explains the eccentricity gap. Note however the sur- 
prising exception of GW Orionis, a system still embedded 
in a disc of Md ~ 0.3 M0, with a period of 2 42 d and an un- 
expec ted low eccentricity of e w 0.04 ±0.06 (jMathieu et al.l 
Il991h . Following the evolution of pre-MS with disc-binary 
interaction is however not the scope of this paper. 

Note that an eccentricity gap seems also present in post- 
AGB binaries (see Fig. [3]) with periods longer than 10"^ d, 
in which tides are inefficient. 



3. Modelling post-AGB binaries with circumbinary 
discs 

We introduce the CB-disc properties and apply our model 
of Section [^U to derive the evolution of post-AGB systems 
in the period-eccentricity plane. 



3.1. Circumbinary disc model 

The post-AGB disc properties can be derived from 
spectral-energy-distribution (SED) modelling in the near- 




Figure 2. Disc geometry as described in iDullemond et all 
(2001) . Li, L2 and L3 are the Lagrangian points, where 
forces cancel out. rin, R and H are the inner radius, the half 
angular-momentum radius and the thickness of the disc, 
respectively. The given radii are indicative. 



Parameter 


Range 


Adopted 


Md/Mo 


10-'' - 10-^ 


10-^ 


td/yr 


2.5 X 10'' - 10^ 


2.5 X 10-^ 


a 


0.01-0.1 


0.1 


H/R 


0.1-0.2 


0.1 


Tin 




Eq.ITO] 


rout/AU 


100-2,000 


500 


/i/Mo 




0.3 



Table 1 . Possible range and adopted values of the param- 
eters, where A/^, td, a, H/R, r^^ and Tout are the mass, 
the lifetime, the viscosity parameter, the thickness, the in- 
ner radius and the outer radius of the circumbinary disc, 
respectively (see also Fig. [2]). /x is the reduced mass of the 
binary. See Section [3. II for more details. 



IR. Observed disc mass es range from about 10 ^ to 10 ^ 
M0 (jGielen et al.l |2003) ■ The rel ative thickness near th e 
inner edge is H/R = 0.1 - 0.25 (Du Uemond et al.ll200lD . 
However, as our model is only valid for thin discs, we as- 
sume H/R = 0.1. 

The inferred location of the inner disc edge mainly de- 
pends on the opacity of the grains which in turn is related 
to their chemistry and size distribution. Metallic iron has 
no spectral feature but a large opacity so, depending on 
the adopted abundance of iron, the SED can be fitted with 
different inner radii ranging from about 2 to 10 AU. These 
values are in good agreement with our predicted inner ra- 
dius (Eq. ITU)) assuming typical separations of post-AGB 
systems from about 0.5 to 5 AU. The out er radius ranges 
from about 100 to a few thousands AU (|de Ruvter et al.l 
2006) and is taken to be rout — 500 AU. The surface den- 
sity distribution decreases with distance fr om the centre of 
the d isc, a{r) ex r^ , with —2 < (5 < — 1 ()de Ruvter et al.l 
12006"). In our model, we choose a{r) ex r~^ such that the 
half-angular-momentum radius is i? = y^Thlrout ■ 

Since when e = the pumping of the eccentricity is 
ineffective (e = 0), we assume that perturbations of the 
orbit, e.g. due to the CB disc or to stellar pulsation lead to 
a minimum eccentricity of 10~^. 



3.2. Post-AGB evolution in the e 
interaction 



logP plane with CB disc 



In Fig. |3| we show four evolutionary tracks in the e — log P 
plane of post-AGB systems with initial periods of 100 d. 
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Figures. Evolutionary tracks during the post-AGB phase 
(black lines) in the e — log P plane including the interac- 
tion with a circumbinary disc as described by Eqs.[7]and[5] 
(see also Fig. [T]) with the adopted circunibinary-disc prop- 
erties of Table \T\ The o symbols are observed post-AGB 
systems. The blue dashed lines show the final eccentricity 
corresponding to eccentricity-pumping as given by Eq. ^ 
modulated by the labelled values (1 corresponds to the 
adopted parameters). The timescale needed to pump the 
eccentricity can be estimated from Eq. [T^l 



300 d, 1, 000 d and 3, 000 d, which interact with a CB disc 
with properties as given in column 3 of Table [T] The binary 
systems are characterised by a reduced mass /i — 0.3 Mq. 
Their orbits are initially circular because tides are very effi- 
cient during the AGB phase and circularise all systems with 
P < 2,0Q0d. Longer-period systems can remain eccentric 
at the start of the post-AGB phase. The blue dashed lines 
indicate the final binary eccentricity due to eccentricity- 
pumping as given by Eqs.[7]and|n]for different disc proper- 
ties. Each line corresponds to pumping modulated by the la- 
belled value. For example, the line labelled 0.5 corresponds 
to the case where e is decreased by a factor 2 compared to 
our adopted case (labelled 1). 

The mass and lifetime of the disc are Md — 10~^ M0 
and td = 2.5 x 10"* yr, respectively, which correspond to 
the highest mass observed and to a lower limit of the disc 
lifetime. As expected from Eqs.[7]and[Sl the e— pumping de- 
creases with increasing period because of the dependence 
of e on the orbital period, i.e. e oc ailt, ex P^^'^. However, 
this trend is not observed among the post-AGB stars (see 
discussion in Section r4.2p . Nevertheless, our simulations can 
reproduce the observed eccentricities of all post-AGB sys- 
tems, except for the three systems with P > 10'^ d and 
e > 0.4. Moreover, all the systems are not equally evolved 
on the post-AGB phase, so some systems may have inter- 
acted longer with their CB disc leading to a more eccentric 
orbit. When extending the disc lifetime to 10^ yr, its ex- 
pected maximum lifetime, our model predicts systems with 
e = 0.7 at P ~ 1Q3 d (see Eq. [T2]). 



4. Discussion 

4.1. Model uncertainties 

The time to pump the binary eccentricity from to e/ is 



t{ef) = 2 X lO^e^. 

fi IO^^Mq 0.1 / 0.1 



0.3 Ma 



H/R 



1000 d 



2/3 



rout 
500 AU 



yr- 
(12) 



given by 



Regarding the possible range of the disc parameters (i.e. 
M(j, td, a, H/R and Tout) given in Table[TJ the eccentricity- 
pumping efficiency can be increased by a factor 36 or de- 
creased by a factor 2000 compared to the adopted values. 
Fig. [3] shows that for an efficency decreased by a factor 
more than 100, the CB disc has a no significant effect on 
the binary eccentricity any more. 

4.2. Discussion 

Our model is able to reproduce most orbits of post-AGB bi- 
naries except for the three systems at P > 10'^ d and e > 0.4 
(Fig. [3]). Extending the disc lifetime to 10^ yr, the expected 
upper limit, gives enough time for the CB disc to increase 
the system eccentricity to the observed value. However, sys- 
tems with periods longer than about 2, 000 d are not ef- 
ficiently circularised during the AGB phase, so they can 
start the post-AGB phase in an eccentric orbit. Moreover, 
as the CB disc is formed during the AGB phase, eccentric- 
ity pumping takes place before the post-AGB phase. This 
leads to an underestimation of the binary eccentricity at 
P > 2,000d. At shorter periods, binary systems enter a 
CE which circularises the orbit before the system starts its 
evolution on the post-AGB phase. 

There is an important discrepancy between the expected 
and observed distributions of post-AGB systems in the 
period-eccentricity plane. Our model predicts systems with 
eccentricities that decreases with increasing periods while 
the observed distribution shows the opposite. As discussed 
before, this discrepancy may partly be explained by an un- 
derestimation in our model of the eccentricity at long or- 
bital periods. A hypothesis to explain this discrepancy is 
that the properties of the disc and the binary are correlated 
in such a way that e-pumping is more efficient in longer- 
period systems. It would be a surprise if discs that form 
through wind mass-transfer or remaining from CE ejection 
are similar. However, our understanding of disc formation 
is poor, so no a priori correlation can be suggested. None of 
the observed disc properties are correlated with the binary 
properties. Nonetheless, it is surprising that no post-AGB 
stars with a disc are observed in systems with orbital peri- 
ods shorter than 100 d. From binary evolutionary models, 
such systems are predicted to exist as post-CE systems and 
are expected to have periods from a few hours to hundreds 
of days. As none are observed, CB discs may be very short- 
lived or unstable at periods shorter than about 100 d. 

In our model we do not consider accretion of matter onto 
the central stars or outflow from the disc. While the effects 
of such re-accretion are observed in some post-AGB stars 
(Waters et al. 1992). our model disc mass is assumed con- 
stant for simplicity. However, as discussed in Section 12.31 
the post-AGB systems are observed at least 10^ yr after 
they left the AGB phase. The inferred disc masses thus 
represent a lower limit of their initial masses. Note that the 
most important effect of re-accretion is to slow down the 
evolution along the post-AGB phase. It then takes longer 
for the central star to become hot enough to efficiently evap- 
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orate the disc gas due to ultraviolet radiation compared to 
single star evolution where no accretion takes place. 

The conclusions reached in this paper not only concern 
post-AGB systems, but a variety of classes of post-mass 
transfer systems such as barium stars, S stars, subgiant 
CH and CEMP binaries, which we observe long after mass 
transfer finished. If CB discs play an important role in the 
evolution of post-AGB systems they are also relevant for 
these classes of stars. Their similar periods and eccentric- 
ities to post-AGB stars support an identical eccentricity- 
pumping mechanism. We are currently working on popu- 
lation synthesis models including the formation and inter- 
action with a CB disc, which attempts at explaining the 
orbital properties of all these systems (Dermine et al. 2012, 
in preparation). 



5. Conclusions 

Our circumbinary-disc model describes the resonant inter- 
action between a disc and its central binary with which we 
account for the large eccentricities observed among post- 
AGB stars. Stable dust discs are detected in a large frac- 
tion of post-AGB stars and are known to be closely re- 
lated to binarity. Circumbinary discs strongly alter the 
evolution of the binary system because of resonant in- 
teractions and reaccretion from the circumbinary disc. 
Resonant interactions with a circumbinary disc transfer an- 
gular momentum from the binary orbit to the disc and 
increases the binary eccentricity . We use the model de- 
scribed bv lLubow fc Artvmowiczl (|1996f ) with the disc prop- 
erties derived from the observed SED of post-AGB stars. 
Reaccretion of gas deficient in refractory elements, because 
they condense in grains maintained in the disc by radiation 
pressure, is responsible for the depletion patterns observed 
at the surface of numerous post-AGB systems and slows 
the evolution of the post-AGB star. 

We estimate binary post-AGB lifetimes to range from 
2.5 X 10"* yr to lO'^ yr which are long compared to single-star 
models that predict only 10^ to 10^ yr. Altough the disc 
lifetime represents the major uncertainty, our model repro- 
duces the eccentricities of most post-AGB systems on the 
derived lower limit of the disc lifetime. Three long-period 
systems show higher eccentricities than expected and re- 
quire more efficient eccentricity-pumping. However, extend- 
ing the disc lifetime to its expected maximum lifetime of 
10^ yr solves the problem. On the other hand, our model 
underestimates the eccentricities of long-period post-AGB 
stars because we assume all systems to have initially circu- 
lar orbits while systems with periods longer than 2,000d 
are not efficiently circularised by tides during the AGB 
phase. Moreover, the circumbinary disc is formed during 
the AGB phase so that eccentricity-pumping already oper- 
ates on the AGB. 

Finally, Ba, S, CH and CEMP stars experience a simi- 
lar evolution than post-AGB stars and show similar period 
and eccentricity distributions. This suggests that interac- 
tion with a circumbinary disc is also a key mechanism in 
the evolution of Ba, S, CH and CEMP stars. 
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